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Bubble collisions in cosmological phase transitions are explored, taking the non-abelian 
Q\ . character of the gauge fields into account. Both the QCD and electroweak phase 

transitions are considered. Numerical solutions of the field equations in several limits 
are presented. 
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Introduction 



The investigations reported in this talk have been motivated by an interest in studying 
cosmological phase transitions quantitatively, taking the non-abelian character of the gauge 



?_i ■ fields into account. Ultimately, we hope to identify observable consequences of cosmological 



phase transitions. 

First-order phase transitions proceed by nucleation of bubbles of the broken phase in the 
background of the symmetric phase. Bubble collisions are of special interest, as they may lead 
to observable effects such as correlations in the cosmic microwave background (CMB)[3] or as 
seeds of galactic and extra-galactic magnetic fields [2] . The quantum chromodynamic (QCD) 
and the electroweak (EW) phase transitions are both candidates of interest in these respects. 
The Lagrangian driving both the QCD and the EW phase transitions are essentially known 
and make it possible to approach the physics of the phase transitions from first principles. 
However, a difficulty to making reliable predictions is that the fundamental guage fields in 
both these instances are non-abelian: the gluon field in QCD and the W and Z fields in 
the EW case. The quantitative role of non-abelian fields in cosmological phase transitions is 
poorly known and difficult to calculate due to the nonlinearities arising from the non-abelian 
character of the gauge fields. 
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2 A Toy Model 



The results reported in this talk are preliminary and correspond to equations of motion that 
follow from toy model Lagrangian, 

£ = £«+£( 2 ) (1) 

£W = — Wt v W iiiV - -B^B^ 
4 ^ 4 M 

9 9 

with 



£ (2) = |(^- 7 |t.^-|s,)$| 2 -f($) 



= w-w;-^#X c (3) 

B^v = d^B v — d v B^ . 

The potential V($) in Eq. leads to spontaneous symmetry breaking that will drive the 
first-order EW phase transition of interest in this work. The detailed form of V(<J>) depends 
upon the theory, but one possible form is given by 

V($) = V(f)= -^ 2 |$| 2 + A|$| 4 . (4) 

This Lagrangian consists of two coupled sectors. One of these is abelian, consisting of the 
scalar boson $ and the vector field B®, which are the familiar scalar Higgs and the vector 
boson of the weak hypercharge current, respectively, of the Weinberg-Salam modelj^]. A first- 
order EW phase transition requires that the Lagrangian be considered in the framework of 
the minimal supersymmetric model (MSSM) extension of the Weinberg-Salam model, and 
one may eventually need to explicitly incorporate the stop field for consistency, as discussed 
in more detail in Ref. pQ. The other sector is non-abelian and consists of the vector field W^. 
We will interpret this field as the gluon field of SU(2) co i or in our study of the QCD phase 
transition in Sect. 3 and as the three vector gauge bosons of the Weinberg-Salam model in 
our study of the of the EW phase transition in Sects. 4 and 5. This interpretation is possible, 
in the spirit of our toy model, since both fields satisfy the same field equations in the absence 
of coupling (7 = 0) between the sectors. Fermions do not appear in the toy Lagrangian, 
which reflects our focus so far on issues of non-abelian dynamics that arise in their absence. 

Equations of motion are obtained by minimizing the action, 

S[ J d A x(C w + £ (2) )] = . (5) 

The result of doing this yields two 'W-equations" : 

d 2 Wl - d,d u Wl - ge 3jk [W*d»Wi + 2W^W* - WfaWft) (6) 

+g 2 e 3jk W^ lm WlW? - f{x) 2 ig 4> v {x) = 0, i = 3, 
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d 2 Wl - d^W; - ge ljk [WXwi + W^Wt - WjftWft) (7) 

+g % k Wie klm WlW™ + lf{x) 2 g 2 Wl = 0, z = (1,2), 



a "0-equation" , 
a "5-equation" , 
and an "/-equation' 1 



d a f{x) 2 ^ a (x) = 0, (8) 
d 2 B u - d^dvB, + f{xfg'^ v {x) = 0, (9) 



^n 2 r)V 
d 2 f(x) - l^f{x)[{Wlf + (W 2 ) 2 ] - f(x)ij a r + f(x)^ = . (10) 



The quantity ip a is denned as 

^ a (x) = d a Q + 9 -B a - ^-Wl . (11) 

3 Bubble Collisions in QCD 

As discussed in Sect. 2, the equation of motion for the gluon field in SU(2) co i or is the 
same as that of above in the absence of coupling between the abelian and non-abelian 
sectors, i.e 7 = 0. In the gauge d^A^ = 0, this is 

d^A a u + ge abc (2A b l d^A c u - A^A^) + gh abc e ceJ ' A b ^ e A{ = 0. (12) 
This equation has a well-known non-perturbative solution, namely the BPST instantonjl], 

a;(x) = (is) 

p g x z + p z 

Using an instanton-like ansatz in 1+1 dimensions, 

A%{x,t) = \ ailv x v F{x,i) , (14) 

bubble collisions were found in Ref. to evolve as solutions of the equation 

d 2 F = —d x F - 12F 2 + 8(x 2 - t 2 )F 3 (15) 
x 

in Minkowski space, assuming periodic boundary conditions in time with 

= 7 — ^T—^+ f r^WT 2 ( 16 ) 

(x — 3) z + p l [x + 3) z + p l 
d t F{x, 0) = . 

The numerical solution to Eq. (fTB^) is given in Fig. ^ Note the development of a gluonic 
wall at the x=0 collision region. The wall grows rapidly beginning at time t = 1.1, and due 
to the singularities in the solution the accuracy of the calculations for t > 1.0 is limited. A 
possible connection to CMB correlations is discussed in Ref. pp. 



3 



Figure 1: F(x,t) for the two bubbleswith instanton-like form 



4 Bubble Collisions in the Abelian Higgs Model 

The abelian Higgs model has been of interest as a prototype for the generation of magnetic 
fields in the early universe in collisions of bubbles during a first-order EW phase transition [2J 
IH1 IE] • The Lagrangian of the abelian Higgs model describes a complex scalar field coupled 
to the electromagnetic (em) field A^ m . It is defined by the Lagrangian Eq. (0) in the the 
abelian sector identifying A e ™ with the field and the electric charge e with the coupling 
parameter g' as e = g'/2. Equations of motion for the Higgs and magnetic field may be read 
off the results in Eqs. (IHl EHTUj) for 7 = 0, i.e., 

d a f(x) 2 M%) = , (17) 
d 2 B v - d»d v B^ + 2e/(x) 2 ^(x) = 0, (18) 



and 



The quantity ip a is now 



dV 

d 2 f(x) - f(x)i> a r + f(x) — = . (19) 



df 



M%) = d«Q + eA e a m . (20) 

In this case, the phase transition is driven by the dynamics of the Higgs field along the 
lines studied by Coleman[7j. Coleman studied the case of a real scalar field (0 = 0), in 
which case the equation of motion for the scalar field becomes 

r)V 

d 2 f(x) + f(x)— = 0. (21) 

The potential V(f) given in Eq. (j3J) describes the dependence of the energy of vacuum on 
the scalar field. It has two minima, corresponding to "true" and "false" vacua. In Coleman's 
model, a symmetry breaking term is added to V(f) to give the true vacuum a slightly lower 
energy. 

One imagines that the system (here, the universe) begins in the false vacuum, and then, 
as time evolves, a transition is made to the true vacuum. The phase transition proceeds as 



bubbles of the true vacuum nucleate in the false vacuum. Nucleated bubbles are tunneling 
(instanton) solutions of the /-equation in Euclidean space[7j. Once nucleated, bubbles grow 
and collide as Minkowski space solutions of the /-equation. 

In their analysis of the abelian Higgs model, Kibble and Vilenkin[2] suggested one way 
in which magnetic fields might be generated as bubbles collide. They considered the regime 
of gentle collisions, where f(x) remains constant, or nearly constant, in the region of overlap 
of the colliding bubbles. They were able to gain insight into the generation of magnetic 
fields in this case by making an expansion about point f(x) = /o, which we shall refer to 
as the Kibble- Vilenkin point. For the case of gentle collisions, we can assume the following 
expansion, 

f(x) = fo + aSf{x) (22) 
if> a (x) = a^ (1) + a 2 V>i 2) + ... , 

where a is the magnitude of fo — f{x) and is small by assumption. Substituting the expansion 
into the equations of motion and requiring that the equations be satisfied at each order in 
the expansion parameter a, the relevant and B equations give, to leading order in a, the 
following results 

(<9 2 + 2e 2 /o 2 )^ 1} = (23) 
<9>« = , (24) 

where now 

= d a Q + eA e a m{1) ■ (25) 

These are essentially the equations of Kibble and Vilenkin who demonstrated from them 
that magnetic fields encircle the overlap region of the colliding bubbles when the phase of 
the Higgs fields is initialy different within each bubble. Their analysis has been elaborated 
upon by Copeland, Saffin, and Tornkvist [H] , who presented solutions in a convenient, closed 
form. 

For violent collisions, where f(x) changes substantially in the collision, the character of 
the problem requires numerical integration of the full coupled PDE. We close this section 
by showing numerical solutions of Eq. (|21jl using an algorithm for solving the PDE in 2 + 1 
dimensions. The initial condition is shown in the left-hand panel of Fig. El and the solution 
of the equations of motion just after the collision begins is shown in the right-hand panel. 
The solutions have been followed sufficiently far in time to convince us that the algorithm 
is stable, and it is easily generalized to include multiple coupled fields. 

5 Bubble Collisions in the EW Phase Transition 

One of our main interests is to understand the generation of the em field during the EW 
phase transition, specifically when non-abelian gauge fields play a role. In this case, 7=1, 
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Figure 2: Two bubbles colliding in the Coleman model as explained in the text 



and the W, B, and <ft (Higgs) fields are fully coupled. The physical Z and A em fields are 
determined in terms of these fields as 

A T = y^yf^'W^gB,) (26) 

^ tf + gB^ W'-gBj. (27) 

We consider here only the case of gentle collisions, extending the analysis given by Kibble 
and Vilenkin for the abelian Higgs model. In the future we will examine the case of more 
violent collisions solving the equations of motion numerically. 

As in the abelian Higgs model, the MSSM equations simplify upon expansion about the 
Kibble- Vilenkin point, 

f(x) = f + a6f(x). (28) 

The fact that ip and W d (for d = (1,2)) enter quadratically in the /-equation places two 
important constraints on these quantities: (1) ip and W d (for d = (1,2)) must have an 
expansion in odd powers of a 1//2 , if we require the square of these quantities be analytic in 
a; and, (2) expanding this equation to leading order in a 1 / 2 , we find that the terms ip^°\ 
w^ 1 , and w^ 2 must vanish. This is most easily seen in the Euclidean metric, from the fact 
that the square of each enters with the same sign. However, the same must be true in the 
Minkowski metric as well by analytic continuation. In view of these considerations, ip a and 
W d for d = (1, 2) have the following expansion 

Ms) = « 1/2 ^ (1) + « 3/2 ^ 3) + (29) 
W d = a^wW + a 3 l 2 wf d + ... . (30) 

It is natural that an expansion in the same parameter a 1 / 2 remains appropriate for d = 3. 
However, there is no requirement that the leading term vanish, so we take 

Wl = w^ 3 + a l ' 2 w^ + a 3 / 2 ^ 3 + ... . (31) 

The B-, 0-, and IV-equations then give, to first order in a 1//2 , 

[d 2 + ^(g 2 + g' 2 )]^ = (32) 
<9>i 1} = (33) 



6 



where now 



1>2Hx) = a a e - £±f^ zv . (34) 

Comparing these equations to those in the abelian Higgs model, Eqs. (|23I24I25|) . we see that 
the Z field here plays the same role as the em field did in that case. Specifically, in this case, 
the phase difference of the Higgs field now determines the gauge field Z M of Eq. (|27jl within 
the bubble overlap region. Thus, for gentle collisions, the mathematical problem in leading 
order is no more complicated that it was in the abelian case. 

We may obtain A em from Maxwell's equation, formed by taking the linear combination 
of the W^- and S-equations, Eqs.fjHJ) and 0, suggested by Eq (j^Bj). We find that the first 
non-vanishing contribution to the em current occurs at order a 3 / 2 and depends cubically 
upon the non-abelian fields Wa^ d calculated at order a 1 / 2 . 

Equations for Wa may be obtained by expanding the B- and jy-equations through 
order a 1 / 2 . We find for d = 3, 

d 2 uf z -d v d-u^ = u (35) 

and 

d 2 u^-d u d-u^ 3 = Sf^ , (36) 

where we have expressed the equations in terms of u v defined as 

u vi x ) = t; w v( x ) • ( 37 ) 

Note that Uv(x) follows from Eq. (j3Tlj) once the driving term ip^(x) has been independently 
determined from Eq. (|34|) and the solution of Eqs. (|32I33J) . For d = 1 or 2 (corresponding to 
d' — 2 or 1, respectively), we obtain the pair of equations 

d 2 u^ d - d v d-u^ d -2[d^u^u^ d ' -u^ d 'uf 3 ) (38) 
+ ( W « d >4°) 3 - u^d^ d ') - (u^ d 'd u u^ 3 - uf 3 d u u^ d ')} 

- 4[(uW 3 ) 2 4 1)d " ^ (0)3 • u^uf 3 } + = . 

Note that Eq. ()38)1 is a linear equation for the non-abelian fields u^' d . The non-abelian field 
u (o)3 en ^ ers nonlinearly, but it is determined from a separate, uncoupled equation, Eq. (|35jl . 
Thus, for sufficiently gentle collisions, all relevant equations are linear. This means, among 
other things, that one can avoid introducing a nonperturbative instanton/sphaleron ansatz 
as in Eq. (|14|) , which would lead to equations that do not transform properly under a lorentz 
transformation. 
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6 Summary and Conclusions 



Methods suitable for exploring the role of non-abelian dynamics in QCD and EW phase 
transitions have been explored. The QCD phase transition is essentially non-abelian, and 
our numerical investigations have been carried out using a BPST instanton-form solution to 
simplify the treatment of the non-linear dynamics. In the case of the EW phase transition, 
we found that for gentle bubble collisions the non-abelian fields may be obtained by solving 
linear equations, so that a nonperturbative sphaleron ansatz is not needed to account for the 
nonperturbative dynamics. Investigations are continuing for the case of gentle collisions along 
the lines outlined here. Our preliminary numerical work for the Coleman model suggests 
that pursuing numerical solutions to the PDE is a promising approach for determining the 
consequences of more violent collisions. 
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